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We first investigate the semileptonic annihilation decay B~ —* rj'l u in QCD. We find 
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the 77' momentum distribution is peaked within its small recoil region due to the loop 
effects. The branching ratio is estimated to be Br(B~ — > rj\~v) = 1.6 x 1CP 4 for 
I = fj,,e, which is accessible at CERN LHC. 
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1 Introduction 



Being consisted of two different heavy flavors, the bottom-charmed meson B c have many 
fascinating properties, which have motivated extensive studies in the literature. Its produc- 
tions HJ, spectroscopy 0,|3|] and decays [|, [5| could be estimated to certain accuracy and 
provid windows for probing both strong and weak interactions. 

The recent observation of B c in 1.8 TeV pp collisions using CDF detector at the Fermilab 
Tevatron has confirmed its existence in nature with mass Mb c = 6.40 ± 0.39 ± 0.13 GeV and 
lifetime t Bc = 0.46ig;^ ± 0.03ps, which agree with the theoretical predictions Further 
detailed experimental studies will be performed at Tevatron Run II and CERN large Hadron 
Collider (LHC). Especially, at LHC with the luminosity C = 10 34 ernes' 1 and ^fs = 14TeV, 
the number of Bf events is expected to be about 10 s ~ 10 10 per year, so that some B c rare 
decays of interests could be studied. 

In this paper, we would like to present the first investigation on the annihilation decays 
B~ — > r]'l~9 in QCD. In this process, b and c annihilate to leptons pair and meanwhile emit 
two gluons to form rj'. Due to OZI suppression, its decay width should be two or three orders 
of magnitude lower than the dominant semileptonic decays induced by b — * c weak current. 
So it belongs to rare decays. However it is still sizable at LHC. 

Compared with the pure leptonic decays B c — > Wi, the suppression factor is a 4 in T(B~ — > 

2 

rjl~v ) instead of the helicity suppression factor — Considering the three bodies phase 
space much smaller than two bodies phase space, one could expect that the decay width 
T(B~ — > rj'l~u ) may be the same order as T(B C — > fiu^) , anyhow it would be much larger 
than T(B C —>■ eu e ). As we know from the radiative J/\l/ decays, the coupling g*g* — > rf is 
much larger then g*g* — > rj, ir°. For example, the ratio BR(J/^> — > r yr]')/BR(J/^> — > 77r°) is 
as large as 10 2 . In the following, the branching ratio is found to be Br(B~ — > rj'l^u) ~ 10~ 4 
for I = fi,e. It is of interests that the common start point suppression factor | k v > | arising 
from the phase space integration d 3 k v > is canceled by the loop functions, so, the distribution 
dBr[B c — > r]'W)/dE v > is peaked within the small recoil region of rj'. This feather makes the 
decay recognizable from the mean experimental background B~ rflu. 
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This paper is organized as followings. In section 2, we give the details of the calculation of 
the amplitude and phrase space integration. Section 3 is devoted to numerical results and 
discussions. 



2 Calculations 

To order of af, the process B~ — > r]'l~u is described by 6 diagrams as in Fig.l. It is 
easy to understood that the heavy quarkonium like B c bound state justifies perturbative 
QCD calculations of the decay. The method outlined years ago for J/\I/ — > 7/7 decays by 
Korner,etaL, |7] is appropriate for the present case. However, we would like to adopt an 
effective Lagrangian approach to avoid introducing B c meson wave function. At first, we 
begin with the sub-amplitude of (be) — > g* a gllv. 



M(bc - glgtl-u) = ^= Vcb g 2 s Tr{T a T b }v c ( Pc 



7a»(1 - 75b j? 7/3-7 77 T° 



+ ^ a u J 1P~P — 2 7n{l ~ 7s) + 7/3—7—77 7m(! - 7s) j 77 7 

fi-Pc-m c fC -p c -m c -p c + fa-m c Pb - fa - m b 



u b (p b ) 

x/y(l - 75)^ + (a «-> 0, h «-> k 2 ) . (1) 



Using the identity 

7 m 7 « 7 /3 = yy*/9 + 7 V° - l a 9^ - ie Ma/35 7575 (2) 

and Dirac equation (p l —m)u(p) = 0, after a bit of algebra, we arrive at an effective Lagrangian 
which has the form 

A = ^V cb g 2 s Tr[T a T b ]c ls (l - ls)b / 7 „(1 - 75)^^ ~ (g^ \ rf). (3) 
Then we can use the definition 

(0\c 1 ,(l- l5 )b\B c (P))=tf B P^ (4) 



and the g*g b — > rf coupling 



isfot I rf) = gXbT^rre^kTk^ (5) 
ki ■ k 2 
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which has been widely used in rf and pseudoscalar productions in heavy quarkonium decays 
and in high energy collidors |J. Here the parameter A v r is understood as a combination 
of SU(3) mixing angles and nonperturbative objects, and can be extracted from the decay 
J/ty — > 7/7. We obtain the total amplitude as 

M = %V A gtTr[TJ' b ]8 ah AA rf if Bc h»(l- lh )p l 



x/2 

1 

2 



x 



2P^ 5 + 2P s K tl - 2P-K 9flS + 2ie wSa P°K< } + AM B m b9llS - Ap blI P s 



K 2 - 2K-p b 



d 4 q 

(2tt)' 



1 



k{ 



I- 6 



1 



k\ — 2p b -k 1 V k\k\ k 2 kyk 2 J k 2 — 2p b ■ k 2 \k 2 kyk 2 k\k 



21,2 



2P^k 5 + 2P s K tl - 2P-Kg^ s + 2ie lipSa P° 'K? + 4M Bc m c g^ s - 4p c ^P s 



K 2 - 2K-p c 



d 4 q 
(27T)* 



kf 
^2u2 



h 2 



k{ 



L.6 * 
~2l2 



k{ — 2p c -k\ \ k(k 2 k 2 kyk 2 J k 2 — 2p c -k 2 \k(kyk 2 k{k[ 



+ 



d A q 2{k 2fl k 2 P ■ k 2 - k lf ,k 2 P-k x 



(6) 



(27r) 4 k 2 k 2 2 kyk 2 \k\-2kyp bi 

where K = k\ + k 2 is the momentum of rf, P is the monetum of B c and q is the loop 
momentum with the relation 2q = ki — k 2 , and the factor | takes into account that both 
sub-amplitudes have already been symmetrized with respect to the two gluons. 

For the heavy b and c quarks, it is reasonable to neglect the relative momentum of the 
quark constituents and their binding energy relative to their masses. In this nonrelativistic 
limit, the constituents are on mass shell and move together with the same velocity. It implies 
the following equations valid to good accuracy 



M(B C ) =m c + m b , p- c 



M 



m b 

Pb = — P ' 

yb M 



(7) 



Hence, we have omitted the terms in eq(6) which are proportional to 



e^Pb^p-c 



It is straightforward to perform the loop integration in eq(6) using Dimensional Regular- 
ization. The amplitude is found to be 

G f 



M = ^V cb gtTr[T a T b }5 ab 4A v af Bc — (P^ + K,f 2 ) l<f(l - l5 )u t 



(9) 



with fi, f 2 defined by 



h = -4C n (K, Pb - K,0,0,m b ) +4C 12 (K, Pb - K,0,0,m b ) 

K K m„' K K m„i 

-2C7ii( —, — - p b , 0, — , m b ) - 2C 12 (— , y - p b , 0, — , m 6 ) 

-4C n (K, Pc - 0, 0, m c ) + 4C 12 (K, p c - K, 0, 0, m c ) 

^ ,K K m n ' \ ^ ,K K mw 

+2C n (y , y - p c , 0, m c ) + 2C 12 (y , - - p c , 0, -j-, m c ) 

2m b K m v > 2m c K m v > 

H Ci2(— ,p& - iv,0, — ,m 6 ) C 12 {—,p c - K,0, —,m c ) 

m c 2 2 m b 2 2 

Ci2(— -p c ,P-K, — Sm c ,m&) 10 

m b m c \ 2 2 

-^ u (y ~p c ,P- A >0,— ,m c ,m b ) 1 , 



and 



— AMrrih ( K K m / 

h = R2 _ 2 K [2C n (K, p b — K, 0, 0, m&) — Ci 2 (K,p b — K, 0, 0, m&) + Cn(—, — — p b , 0, —^-,m b ) 

AMm ( K K m i \ 

+ K 2_ 2 p c .K [ 2Cu Pc ~ K, 0, 0, m c ) - Ci 2 (K, p c - K, 0, 0, m c ) + C u (y , y - p c , 0, m c ) J 

+— fc u (^,P6-^^,0,m 6 ) -2C la (^,»-X > ^,0, m6 ) + <7 (* ^,0,m 6 )) 
m c \ 2 2 22 22/ 

/" 771 f Tfl i Tfl i 

Cn(-,p c -iT,^f Am c )-2C 12 (-, Pc -K,^- ,0,m e ) + C (-, Pe -K,-2- ,0,m c ] 

m b V 2 2 2 2 2 2 

M(m b -m c ) ( K m v > K vn^i 

Cn(- -p c ,P - K, —^,m c ,m b ) - 2Ci 2 (— -p c ,P- K, — <-,m c ,m b ) 



m b m c V 2 ™ ' 2 ' u ' u/ ±iV 2 ' 2 

+c (^-; Pc ,P-K,^,m c ,m b )') (11) 



The scalar loop functions and their definitions can be found in ref ||. The divergences are 
canceled as they should be. 
With eq(9), we get 

dBr(B~ -> Tj'l-v) = 1 f^j — 1(i;T ,, , ...2 , f ,2 r <2. 

dE v > (2tt) 5 1QM Bc 

where 



AirJE*, - m%—Ml c {Eti ~ m\) \fi + f 2 ? C 2 r Bc , (12) 



C = ^a 2 J B A v ^V cb . (13) 
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3 Numerical Results and Discussions 



For numerical results, we would take a s = a s (Ms c ) = 0.2, = 0.04, A v > = 0.2 and 
t Bc = 0A6ps. The decay constant fs c probes the strong(nonpertubative) QCD dynamics 
which bind b and c quarks to form the bound state B c . It is common wisdom to realize 
that the size of B c would much larger than the size of Bq. The size of B c scales as l/m c , 
but the size of Bq scales as l/rriq (q = u, d, s). The compact size of B c would enhance the 
importance of its annihilation decays and imply the decay constant fs c would much larger 
than /b |. In nonrelativistic limit, fs c can be related to the value of its wave function at 
origin JlO] . Using the nonrelativistic potemtal models, Echiten and Quigg |§ estimated 

500 MeV (Buchmiiller-Tye potential [11]) 

512 MeV (power law potential |12|| ) 

479 MeV (logarithmic potential |L3|) 

687 MeV (cornell potential [0]) 

For numerical illustrations, we would take fg c = 500 MeV. The r/ momentum distribution 
is displayed in Fig. 2. We find the rj' momentum distribution is peaked within small recoiling 
region of rf. One can expand the scalar functions in fi and f% in terms of basic scalar 



(14) 



functions Bq and C in [§] and find the factor y-Ey — m^,(E^, —m^t) will be canceled by the 
loop function. For an example, we expand Cn(K,p b — K, 0, 0, m b ) as 

1 



C n (K,p b -K,0,0,m b ) 



x 



2{K\ Pb -Kf-{K-{ Pb -K)f) 
(p b - K) 2 (B (p b , 0, m b ) - B (p b - K, 0, 0) - K 2 C (K,p b - K, 0, 0, m b )) + 

1 r , 



(15) 



x 



Therefore, the rj' momentum distribution would behave as 

1 

oc 



(16) 



when E v > is small. The singularity at the start point of the distribution due to the factor in 
eq(15) is integratable and give finite decay width. Such peculiar property would make the 
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decay itself recognizable from its mean background B~ — > rflu at LHC, especially, when the 
decay chain rf — > 77 could be used to reconstruct the events in data analysis. 
The branching ratio is estimated to be 

Br{B~ -> rfrv) = 1.6 x 10" 4 , (17) 

which is accessible at LHC. We can extend the estimation to Br(B c — > 7T W), if the following 
relation is valid 

Br(B c -> rfiv) Br(J/^ -> 7/7) „ , s 
L_£ ' — L Li L^~10 2 , (18) 



5r(5 c -> tt°£p) Br(J/^> -> 7r 7) 
which is the ratio square of the relative 7/ and 7r° coupling strengths to gluons, when the 
phase space diference is small due to m v ' and m n o. In this way, we get 

Br(B c -> 7T°£z/) ~ 1.5 x 1(T 6 . (19) 

In conclusion, we have presented the first study on the semileptonic annihilation decays 
B~ — > v(l~v . The rj' momentum distribution in the decay is found peaking within the small 
recoil region of 77' because of the loop effects, which is different from the common tree level 
cases. The branching ratio is estimated to be few times larger than the pure leptonic decays 
B~ — ► \w and 10 4 times larger than B~ — > ev. With large samples to be obtained at 
LHC, the decay could be measured and might be used to extract B c decay constant and/or 
to probe strong and weak interactions. 
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Figure Captions 



Figure 1: Diagrams for B c — > rj'l v at the leading level. The blob represents r/' 



Figure 2: The distribution of dBr(B c — > rj'l u)/dE n i as a function of E v / 
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